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A heuristic theory for VLF and LF fields near and in a rough sea surface is obtained by 
first finding the electromagnetic field configuration in air and then calculating the undersea 
fields by means of the Helmholtz-KirchofF Integral Theorem. The fields above the water 
are found by a succession of quasi-static approximations which depend on the observation 
that the scale of irregularities on the sea surface is very small compared with an 
EM wavelength. _> 

The theory predicts that the configuration of the //-field above the water depends on 
the direction of EM propagation relative to the wave crests. It also predicts that for 
underwater measurements made a few tens of feet below the troughs of the waves, the field 
variations due to one- or two-foot sea waves are averaged out; but for storm waves the 
phase and attenuation of the field observed underwater varies with instantaneous water 
height. These theoretical predictions have been confirmed experimentally. 



1. Introduction 

The well-known stability of LF and VLF propaga- 
tion at great distances has produced an increasing 
number of applications to international comparison 
of frequency standards and to navigational aids. 
In this latter connection, it is important to have 
some estimate of local perturbations in field caused 
by a rough or stormy sea surface. Since at these 
frequencies, the attenuation of radio waves is only 
of the order of a few decibels per foot in sea water, 
it is also relevant to try to estimate the fields meas- 
ured by an antenna located a few feet under the 
water surface. 

Both the above-the-surface and subsurface aspects 
of this problem have received some attention [Rice, 
1958; Herman, 1958; Wait, 1959] in the last few 
years. Only one of the works referred to, that of 
Wait, treats the problem of subsurface fields. Wait 
uses the "Leontovitch" [1948] assumptions concern- 
ing the relative magnitudes of the radius of curvature 
of the sea surface, the inverse of the propagation 
constant in sea water and the distance under the 
sea surface at which one wishes to calculate the 
fields. 

There are, however, naturally occurring situations 
in which the Leontovitch assumptions do not hold, 
and in which we still may wish to calculate the fields. 
Our approach is to find a solution to the field prob- 
lem above the sea surface on the assumption that 
the sea surface is perfectly conducting (to the extent 
that the E-field is everywhere normal to the surface) 
and that the water waves make the sea surface a 



i Operated with support from the U.S. Air Force. 



trochoidal cylinder; then (dropping the assumption 
that sea water is a perfect conductor) we use the 
integral theorem of Helmholtz and Kirchoff to find 
the fields under the sea surface. 

The general behavior of the fields above the sea 
surface can be readily ascertained from the equations 
which are developed for these fields. The E-field 
in the air is locally enhanced over the crests of the 
waves and diminished over the troughs, in accordance 
with the local curvature of the surface. Whether 
this effect is passed on to the //-field and (through 
continuity of II at the boundary) to the subsurface 
fields, depends on the angle between the direction 
of propagation and the direction of wave crests. 

The behavior of the fields at some distance 
beneath the sea surface is not so evident from the 
equations derived for them. The theoretical re- 
sults obtained have been evaluated for typical sea 
conditions at VLF. The result of these calculations 
is to predict that an observation made at reasonable 
depths under storm seas should show variations in 
amplitude and phase proportional to the instan- 
taneous height of water above the point of observa- 
tion. On the other hand, in calm seas (waves up 
to 1 or 2 ft) the effect of the waves on the field ob- 
served at depth is averaged out. 

2. Trochoidal Coordinates 

The shape of a periodic gravity wave on an in- 
finitely extensive sea tends to be a trochoid [Lamb, 
1945], a curve which is a generalization of a cycloid. 
Such curves are described by a point on the radius 
of a circle as the circle is rolled on a straight line; 
they tend to have steeper peaks and broader troughs 
than would be the case for a sine wave. 
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Actual sea waves depart from the cycloidal curves 
assumed in this paper in two ways: First, if the 
waves have sufficient amplitude so that they are 
cusped and breaking, the cycloidal-type trochoid 
assumed is only an approximation to the actual 
surface [Lamb, 1945]. Second, in actuality the 
wind excites a random spectrum of waves, so that 
each "wave" is actually a superposition of wavelets 
and no two waves are really alike. Moreover, the sea 
is sometimes "crossed" by waves propagating in 
several directions. Further, although real sea wave 
"crests" do stretch for distances comparable to a 
wavelength, they do not stretch indefinitely to pro- 
duce a cylindric sea. 

The willingness to ignore these differences be- 
tween the real sea and the model is based on the 
assumption that local sea-surface conditions affect 
the surface fields in that same local area only. 
For the fields above the water, this assumption 
can be justified on the grounds that local surface 
irregularities (the water waves) would produce 
nonlocalized effects only if they tend to scatter the 
propagating radio wave. These irregularities, how- 
ever, have dimensions which are so small compared 
with a radio wavelength in air 2 that little such 
scattering can occur. 

Beneath the sea, the assumption of local influence 
is justified by the high attenuation of radio waves 
propagating through the water. The field beneath 
the sea surface can be thought of as due to a con- 
tinuous distribution of sources along the surface 
(the strength of the field to be calculated by a 
Green's Theorem approach or the Helmholtz- 
Kirchoff Integral Theorem). The effect of a given 
source area diminishes exponentially with the distance 
between the surface source and a given subsurface 
observation point. Since this attenuation is of the 
order of 1 dB per foot at 15 kc/s, the subsurface 
field is for all practical purposes determined by the 
surface field in the local area above the subsurface 
point. 

Thus, if we find the infinite trochoidal cylinder 
which matches the local sea conditions we are 
concerned with, we may trust that the failure 
of the match in other areas of the sea surface will 
not produce local fields noticeably different from 
these which we calculate. 

Since the boundary conditions we must satisfy 
are applied on the surface of a trochoidal cylinder, 
it is convenient to use trochoidal coordinates in 
solving the field problem. We generate these 
coordinates as follows: 

A A 

let z=x+jy w=u+jv, j=\ — l (2.1) 

A A 

z and w represent points in two Cartesian planes, 
the abscissa and ordinate being given, respectively, 

A A 

by x and y in the z plane and u and v in the w plane. 



- At 15 kc/s the electromagnetic wavelength in air is 20 km. Ocean wave- 
lengths seldom exceed 300 m. Ocean wave heights (peak to trough) seldom 
exceed 10 m. 



Further, let 

A A A 

w=z+ja exp(jz). (2.2) 

Equation (2.2) defines a transformation from the 

A A 

z plane to the w plane. All straight lines parallel 

A 

to the x axis in the z plane are transformed or 

A 

mapped by (2.2) into trochoids in the w plane 
(whether or not ae~ v <i\). The mapping of the 

A A A 

z and w is also conformal for values of z such that 
ae~ v is less than unity. We now define trochoidal 
coordinates 6' and n, in the w plane, as follows: 



n = y. 



(2.3) 



The curves, n = const are trochoids, and further- 
more if i/>ln a, 3 the curves s= const ant intersect 
the curves n= constant at right angles (because the 
transformation (2.2) is conformal). 

We may combine (2.1) and (2.2) to obtain 



u=x—ae~ y sin x 
v—y-\-ae~ u cos x 
and using (2.3) we get 

u=s—ae~ n sin s 
v=n-\-ae~ n cos s. 



(2.4) 



(2.5) 



The scale factor h (the same for n as for s because of 
conformality) may be found from 



so that 



h=(l—2ae- n cos s+a 2 e~ 2n ) 



1/2 



(2.6) 
(2.7) 



We need only add a third coordinate axis normal 
to the u, v plane to change trochoidal coordinates 
to trochoidal cylindrical coordinates. If the third 
coordinate is denoted by z, then u, v, and z are 
the Cartesian coordinates of the three-dimensional 
space. The scale factor associated with the third 
coordinate is, of course, unity. We will denote 
the unit vectors in the s, n, and ^-directions by 

e s , e n , and e z . 

3. First Solution in Air 

A radio wave propagating over the surface of an 
unrippled sea at or below broadcast frequencies is 
for all practical purposes vertically polarized. To 
find the field outside the sea, we assume the sea to be 
a perfect conductor with the surface having the form 
of a trochoidal cylinder. 



3 We use In to denote the natural logarithm. 
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The conductivity of the sea water is so large (4 
mho-meter" 1 at standard temperature and salinity) 
compared with the susceptance (1()~ 5 mho-meter" 1 
at about 150 kc/s) of the air above it, that what little 
loss power that enters the sen propagates essentially 
straight down. If the sen surface is disturbed by 
waves, the EM field above it will be correspondingly 
perturbed. Since the conductivity of even a few 
inches depth of water is high compared with either 
the susceptance or the Y of the air immediately 
above it, we can expect that the effect of the disturbed 
sea surface on the E'-fields above it will be nearly the 
same as that of a perfect conductor — the externa] 
^-fields must remain perpendicular to the sea surface. 

In the case of a static I^-field in cylindric coor- 
dinates, the problem is one of finding a scalar poten- 
tial ip which satisfies Laplace's equation 

vv=o 

and for which <p is constant on the boundary. 

In coordinates derived from a conform al trans- 
formation, the two-dimensional Laplacian V 2 is 
[Morse and Feshbach, 195;]] 



~h 2 W dn 2 / 



(3.1) 



Obviously, then, if <p=E n, <p satisfies Laplace's equa- 
tion and 



now state that 



E{s,n) =—V<p= — -^ e n 



(3.2) 



is the corresponding electric field. 

The equation we are really interested in solving 
is the electromagnetic wave equation. When the 
time variation is separated out of this partial dif- 
ferential equation, the result is called the Helmholtz 

-> 
equation. If F is a vector function of the spacial 
coordinates, the Helmholtz equation has the form 



V 2 F+k 2 F=0 



(3-3) 



where for good dielectrics k 2 =co 2 fji€, and for good 
conductors k 2 =—jw/jLa. Here w is the radian fre- 
quency of the signal; /x, e and a are the (magnetic) 
permeability, permittivity (electric), and conduc- 
tivity of the relevant medium, respectively. 

It is a well-known property [Ramo and Whinnery, 

-> 
1953] of cylindrical systems that if E(£ lt £ 2 ) is a 
static two-dimensional field solution (i.e., one de- 
rivable as the gradient of a potential satisfying 
Laplace's equation) in the plane perpendicular to 
the cylinder axis (£j and £ 2 being the coordinates of 
this plane), then 

satisfies the Helmholtz equation (3.3). One can 



E= 



-E, 



I, 



e- Jkz e n 



(3.4) 



satisfies (3.3) from (3.1), (3.2) and the discussion 
between them. This field expression represents a 
TEM wave traveling along the cylinder axis. In 
the terms of the original physical problem, it repre- 
sents a TEM wave traveling parallel to the sen wave 
crests. If we examine the form of h [see (2.7)], 
we note that the field becomes uniform as n becomes 
infinite. This is what we would expect from physi- 
cal reasoning, as water wave action at the surface 
of the sea should not affect the fields at arbitrarily 
large distances from it. We also note that the rate of 
decay of this "surface effect" has nothing to do with 
the electromagnetic wave number, k. The way in 
which (3.4) has derived from a static field solution 
makes it obvious why this should be so. 
We also note that 



hence 



4. Second Solution in Air 



(3.5) 



The results of section 3 apply only to TEM waves, 
those propagating parallel to t lie wave crests. We 

wish now to examine cases of EM propagation a1 an 
angle with respect to the sea wave crests; in par- 
ticular, propagation perpendicular to the crests. 
This problem is inherently more difficult than the 
previous one, in that it requires manipulation of 
solutions of the vector Helmholtz equation in curvi- 
linear coordinates. Helmholtz' equation for a vector 
-> 
E and a scalar <p are symbolically the same. 



vV+*V=o 



V 2 E+k 2 E=0 



(4.1) 
(4.2) 



but the interpretations that must be given to the 
operators V 2 are (in curvilinear coordinates) different ; 4 
for this reason (4.2) may be solved by separation 
techniques in only six coordinate systems, whereas 
(4.1) can be so solved in eleven. With respect to 
other coordinate systems (including those appro- 
priate to the present problem), the difficulty is not 
so much in finding solutions to the equations, as in 
finding solutions that will also satisfy the boundary 
conditions. 

It is much easier to find solutions to (4.1) than to 
(4.2), whether or not the equations separate. If <p 
is a solution to the scalar wave equation (4.1) in a 
cylindrical coordinate system, then a set of orthog- 



4 vV is div grad ^ for <p a scalar; V 2 E is grad div E— curl curl EforEa vector 
[Moon and Spencer, 1961). 
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onal solutions to (4.2), L, M, and N, may be 
constructed from <p as follows : 



£=: grade 



M==curl (<pe z ) 



iV=curl curl (<pe z ) 



(4.3) 



(4.4) 



(4.5) 



in which e z , the axial unit vector, is parallel to the 

-> 
surface of the water, L is a "longitudinal" wave, 

-> -> 

and M and N are called "transverse" waves. In 

-> 
particular, the tangential component of M can be 

-> 
made to vanish along a boundary parallel to e z ) it is 

therefore the appropriate wave by which to describe 

-» 
the E-field outside the water. Thus, the problem 

-> 
of finding the vector field E reduces to that of finding 
an appropriate scalar <p in the trochoidal curvi- 
linear coordinate system. 

In trochoidal cylinder coordinates, the scalar 
Helmholtz equation is 



^te + » + ^ + ^ =0 - (4 - 6) 



// 



If the propagation is parallel to the wave crests, 
then for plane waves 

dz 2 



and (4.6) reduces to 



ds 2 ^ 



dn 2 



(4.7) 



which is precisely the problem treated in section 3. 
Let the propagation be at some angle 6 with respect 
to the ^-direction; because the boundary has cylin- 
dric symmetry, because the electrical wavelengths 
are very long, and because we would like a solution 
closely related to the plane wave solution of section 3, 
we assume that the variations of <p with z are ex- 
ponential only: 



<p(x, y, z) = {x, y)e~ jkz cosi 



(4.8) 



In this manner (4.6) becomes the two-dimensional 
scalar wave equation : 

KB+») +jvsin ' 9 = o - <4 - 9 > 



In discussing the solution of this equation, nothing 
is lost by setting sin 0=1, i.e., considering propaga- 
tion perpendicular to the wave crests. 

Unfortunately (4.9) does not separate in the tro- 
choidal coordinate system. 5 Even though a suitable 
exact solution cannot be found for (4.9) in closed 
form, it may be possible to build up a series solution 
in ascending powers of k, if k is small ; we can termi- 
nate the series at whatever term yields a satisfactory 
approximation to the desired solution. The method 
of obtaining the series solution is described in 
appendix A, together with the details of the calcula- 
tion for the case of interest here. The result is 



HM 



ka 2 e~ 



-k 2 a(n+\)e~ n 



cos s — 



k 2 a 2 



-**"|- (4.10) 



The principal part of \p in (4.10) is the quasi-static 
term (E /jk) e~ jks . All the other terms are smaller 
by at least a factor of k and die out with increasing 
curvilinear height, n, above the water surface. We 
have assumed a wavelength of 2t for the water 
wavelengths, whose actual wavelengths L range 
between 40 and 300 m. The actual wavelengths X 
of EM waves in air at VLF and LF range from 3000 
to 30,000 m; the corresponding value of k to be used 
in these equations is found by expressing X in units 
of L/2tt: 

thus k ranges from less than 0.001 for short L at 
10 kc/s to 0.05 for long L at 100 kc/s. 

The Afield follows directly from (4.10). 

E 



E 



-VX(tez)=-£e 



..[ 



l-\-jk{a(n-\-l)e w sins}— k 2 



{ 



-e kn +a(n J r l)e n cos s 



yy. 



-JkE , 



h 



jks 



fa 2 . 



-e kn )-\-ane n cos s 



-p.. 

(4.12) 



The main features of this expression are first, that 

the E'-field is essentially a propagating static field; 

-> 
second, that the tangential component of E, the 

-> 
terms in e s , does in fact vanish at the water's 
surface, n=0; third, that the correction terms are 
small. _> 

The expression for the magnetic intensity H is 
somewhat simpler in that we have from Maxwell's 
equations 



H= 



_J_ 

CO/jL 



VXE. 



(4.13) 



5 The condition for separation is that h 2 be the sum of a function of s only and a 
function of n only [see Morse and Feshbach, op. cit. pp. 498-500]. 
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But since div (\pe z ) is zero, we have from Maxwell's 
equations 6 



Hence 



VX#=VXVX(iK)== A 2 (*«,). 



H= J — rpe 2 =jkY xl/e 2 



in which F , the characteristic admittance of free 

space has been written for k/oo/uL. 7 What is note- 

-> 
worthy about this expression for H is that it does 

not contain the scale factor h. Thus, the E fields 
are essentially the same and essentially given by the 
static field configuration for propagation both 

parallel and perpendicular to the wave crests, pro- 

-> 
vided k is small. The 77 fields, however, are not 

-» 
the same; in the former case, the magnitude of II 
varies with the scale factor h from point to point 

near the waves; in the latter case, the magnitude of 
-> 

H is substantially unvarying with respect to position 
above the waves. 

Finally, we return to the question of propagation 
at an arbitrary angle, 0, with respect to the wave 
crests. With respect to the entire argument in 
appendix A, we can replace k by k sin 6 without 
changing either the argument or the result. In 
accordance with (4.8), we can multiply \p(k sin 6) 
by exp (—jkz cos 0) to obtain a function <p which 
satisfies the three-dimensional scalar Helmholtz 
equation 

<p=t(ksmd)e- jkzcosd . (4.14) 



— * 
The S-field is then given by 

E=VXM(k sin S) e~ jkz eoa °e J, 



(4.15) 



but the term exp (—jkz cos 6) is constant with 
respect to the curl operation on a vector in the 
^-direction ; hence 



E= e - jkzcosd VXbP(ksm 6)e z \. 



(4.16) 



Thus the fi'-field is exactly the same as that obtained 
with sin 0=1, except that k is replaced by k sin 6 

and the whole field is multiplied by exp (—jkz cos 0). 

-> -> 

We find II by taking the curl of E to obtain 



H=^-VXE 

LO/JL 



CO/Z 



VX{e- jkzc ° 8d VX[\ls(ksmd)e 2 }}. (4.17) 



e Morse and Feshbach, oj>. cit., p. 1766, eq 13.1.6. 

7 While k or X is measured in units of L/2ir as in (3.41), Y Q is invariant to a change 
of scale in length and is measured in whatever units are consistent with those 
assigned to E and H, e.g., 1/377 mho for E in volts per meter and Hin amperes 
per meter. The units of n are then dependent on those chosen for E, H, and 
length. 



By directly performing the curl operations [Morse 
and Feshbach, 1953, p. 115)] and by noting (4.9), 
we obtain after some manipulations 



H_ _y cos 6 

E h 



9 F dxP(k sin 6) -* b$(k sin 0) ■ 



'«»] 



+Y jk sin 2 d e - Jkzcoad ils(k sin S)e z (4.18) 

in which we have once more written Y Q for &/co/x. 

If we discard all terms having higher order in /*, 

-> -> 

the expressions for \f/, E, and H reduce, respectively, to 

p 1 p — Jk(s sin 9+z cos 0) 



jk sin 6 



Eo 
h 



J£ ^0 e ~Ms sin d+z cos 0) , 



H=Y E (^Z+sm6^ 



(4.19) 
(4.20) 

? -jk(s sin d+z cos 0) (4 21) 



5. Modified Helmholtz-Kirchoff Integral 
Theorem 

The Helmholtz-Kirchoff Integral Theorem [Born 
and Wolf, 1959] provides a method of calculating 
the fields inside a source free volume when the fields 
and their normal derivatives are known on the sur- 
faces bounding it. If a scalar field, U, satisfies the 
Helmholtz equation 



(v 2 +W=o 



(5.1) 



everywhere within a volume V whose bounding sur- 
face is S, and U is continuous and has continuous 
first and second derivatives in V and on S, then the 
field at a point P within the volume is given by 
[Born and Wolf, 1959] 



^JJK(?)-t£> 



(5.2) 



and -r— are taken 
on 



where J is the distance from the point P within the 
volume to the element of area dS on the surface. 

The normal derivatives r— ( — 1 

along the inward normal at the location of the sur- 
face element dS. Although (5.2) is usually derived 
for a scalar field, in Cartesian coordinates, u, v, z, the 

Helmholtz equation for a vector field H in a source 
free isotropic medium is simply 

(^+S+S+ F ) r/ =°- (5 - 3) 

Hence, (5.1) holds individually for each of the 

-> 
Cartesian components of H and (5.2) may be applied 
to each of these in turn. 
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Now, the geometry of our problem is cylindric so 
we may adapt (5.2) to simplify computation in this 
case. We let S be a cylindric surface with the 
cylinder axis in the ^-direction. The Cartesian 
coordinates of the plane perpendicular to the cylinder 
axis at 2=0 will be denoted by u (abscissa) and v 
(ordinate). Furthermore, let us suppose that U(,u, 
v, z) can be expressed as follows: 



U(u, v, z)=U(u, v)e- jbz 



(5.4) 



where b is a constant. Then we may rewrite the 
integral (5.2) as 



uiP) -K[ U{u ' v)e ' 3b t(T 



e jk t dU(u,v) 
£ dn 



-e 



] 



dS. (5.5) 



Now, a plane, z= constant, intersects S in a curve 
C, which is the same for all such constants because 
of the cylindric geometry. In appendix B this sur- 
face integral is reduced to a line integral along C by 
integration with respect to z. The result obtained 
is 



pcn-j/.[i 



THi l) (kr)~-m i) (kr)^ \ 



3< 



(5.6) 



in which r is the vector distance from a point p on 
C to P, C being taken so as to lie in the same plane 

as P. The vector n in (5.6) is a unit inward normal 
from C at point P; the functions H^ and H^ are 
Hankel functions of the first kind, zero and first 
order. This result is the Modified Helmholtz- 
Kirchoff Integral Theorem (MHKIT) for cylindric 
geometry. 

Both (5.6) and the usual form of the Theorem (5.5) 
determine the field U at any point P inside the curve 
C from the value of U and its normal derivative on 
the curve, given k for the material within this 
boundary. While this is indeed a correct result, it 
seems to require too much information in that either 
the field or its normal derivative at the surface 
should suffice to specify the function inside the sea 
[Stratton, 1941]. It is beyond the scope of this 
paper to enter an extended discussion of the mathe- 
matical subtleties of this apparent overspecification 
of the boundary conditions. The following obser- 
vations are nevertheless helpful : 

(a) The Helmholtz-Kirchoff Theorem is essentially 
an application of Green's functions. If we could 
find the Green's function for the bounded region, 
Vj we could calculate the fields in the interior from 
their values on the surface alone [Morse and Fesh- 
bach, 1953, pp. 803-807]. In our problem, we know 
of no way of determining the requisite Green's 
function. The Helmholtz-Kirchoff Integral Theorem 
is a means of using the Green's function appropriate 
to a known solvable problem in a problem which has 
an arbitrary boundary. The usual HKIT makes 



use of the Green's function for a spherical boundary; 
the modified theorem uses the Green's function for a 
circular cylinder. The price that one pays for 
using the simple Green's function is that one must 
know both the field and its normal derivative at the 
boundary. So, if we are able to provide extra infor- 
mation, we can use the MHKIT to avoid mathe- 
matical difficulties. 

(b) Once the field U is specified at the boundary, 
the dUldn required in the equation is also deter- 
mined. It must be calculated or estimated by other 
means before the Integral Theorems can be used. 

The required -^— generally cannot be found from the 

theorem itself, even by successive approximations 
[Franz, 1949; Shelkunoff, 1951]. 

We have already calculated the electromagnetic 
field propagating above a troehoidal sea. 

What remains is the question of determining the 

normal derivative of the tangential H on the sea 
surface. The tangential component of the magnetic 
field is always continuous across a boundary between 
two media, but its normal derivative is not. For 
this normal derivative to be continuous across a 
boundary, it is necessary that both the conductivity 
and permittivity of the two media be the same. 
This condition is not satisfied here. Our method of 
attack is outlined in the next section. 

6. Estimating the Normal Derivative 

Suppose we have a circular cylindrical (air filled) 
cavity in an infinite conducting medium. The con- 

A 

ducting medium has wave number k=(—l+j) 
■y/irfii(T. Now suppose there is a cylindrical scalar 
wave propagating outward from the cylinder axis 
which is given in the conducting medium only by 



U(r)=AHi 1) (kr) 



(6.1) 



where r and 6 are the usual cylindrical coordinates. 
Now, if the radius of the circular cross section of 
the cavity is p, the field on the surface of the con- 
ducting medium may be denoted by U s and we 
have 






(6.2) 



which reduces to U(r) = U s at the cylinder surface 
r = p. For the normal derivative we then have 



dU 
dr 



m»(k P ) 



(6.3) 



What we shall attempt to do is to estimate the 
value of the normal derivative on the sea surface 
by representing the tangential component of the 
field at each point on the surface (or just beneath 
it) by a cylindrical wave of the form given in (6.2), 
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where p will be chosen to be the radius of curvature 
of the cross section of the sea surface normal to the 
wave crests, at the point in question. Hence, (5.6) 
becomes 



U(py- 



kj 
4 



M 



r (kr) '"'' + 1 l£\kr) 



Hl l) {kp) 



1 da. 



(6.4) 



This procedure suggests that the sea surface is made 
to conform locally, with the boundary of the cylin- 
drical cavity referred to above. Now, the center of 
the circle (of radius equal to the radius of curvature 
of the curve at point p) that best conforms to the 
<mrve at p is on the concave side of the curve. 
Throughout most of its length, the cross section of 
the sea surface is concave upward, so the idea of a 
cylindrical wave propagating into the sea from an 
axis above it makes some sort of physical sense. 
On the portions of sea surface cross section which 
are concave downward, we expect the cylinder axis 
in question to be beneath the sea surface. The 
wave appropriate to the physical conditions of the 
problem is then an H ( Q 2) wave. Hi 2) is a Hankel func- 
tion of zero order and second kind, and represents a 
wave propagating inward toward the cylinder axis. 
Now, if 



U(r)=Am 2) (kr)e, 

dU A -* 

dn 



(6.5) 
(6.6) 



The expression has a positive sign because the inward 
normal makes an angle of ir radians with the direc- 
tion of increasing r. If again we are given that the 
surface value of U is U., 



A-- 



U s 



and 



dU 

dn 



H^(k P ) 



' Hi?\kp) 



(6.7) 



(6.8) 



If the rectangular coordinates (u, v) of a curve are 
given parametrically in terms of an auxiliary variable 
x, then the radius of curvature, p, is given by 



m<%)T 



( 



du d?v_ 
dx dx 2 



dv^d 2 u 
dx dx' 



) 



(6.9) 



where the sign is chosen to achieve a positive result. 
In particular, for a trochoid, we can place n=0 in 
(2.3) and (2.5) to obtain from (6.9) the radius of 



curvature 



! 



(1 — 2a cos x+a 2 )* 12 
a(a— cos x) 



(6.10) 



Any point p on the curve, at which the curve is 
concave away from the sea, corresponds to a value 
of x for which (a— cos j)>0. For such values of p, 
we use (6.3) to estimate the normal derivative 
dU/dn. Infinite radius of curvature occurs when 
(a— cos x)=0; and any point p at which the curve 
is concave into the sea corresponds to a value of x 
for which (a— cos x)<^0, for which (6.8) is used to 
estimate dU/dn. 

7. Calculating the Subsurface Fields 

We must now assemble the results derived in 
previous sections to obtain an expression for the 
subsurface magnetic field. 

At VLF and LF frequencies, it suffices to use the 

— » 
quasi-static approximation to H above the water, 
as given by (4.21) 



H=Y E d 



cos 6 



«.) 



e s +sin 6e z )e-J k(8 shie + z G0 



(4.21) 



In order to use the MHKIT on this vector H, we 
need to express it in rectangular coordinates. For 

a trochoid expressed parametrically as in (2.5) with 

-» 
n=0, the unit tangent e s to the water surface is 
given by 

-> 1— a cos x -» a sin x -> , n 1 N 

e s = — t— - e u r — e v . (/.l; 

Thus, the magnetic field becomes in rectangular 
coordinates 



II 



=Y E ( 



(1— a cos x) cos 6 -> a sin x cos 6-* 



h 2 



K 1 



e v 



+ sin o^\-M sshie + zcosd \ (7.2) 

Moreover, s can be expressed parametrically in terms 
of x, so that 

ds—hdx. (7.3) 

Since the MHKIT applies individually to each of 

-> 
the rectangular components of H, it applies to the 

-> 
vector H when it is expressed in rectangular co- 
ordinates as in (7.2). By placing (7.2) and (7.3), 
(6.10), and either (6.3) or (6.8), as required, in (5.6), 
we could write down a general integral for the sub- 
surface field in which all variables are expressed in 
terms of the parameter x. Such a general expression 
is complicated. Either it must simplify by way of 
approximations, or the entire process must be re- 
duced to a computer calculation. Unfortunately, 
at VLF and LF frequencies, skin depths in sea 
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water are of the order of a few meters (1 to 3 dB of 
attenuation per foot in the 10 to 100 kc/s region), 
wave heights are of the order of a few meters, and 
the depth at which one can make observations is 
limited by attenuation from a few meters to a few 
tens of meters. Thus, many of the variables that 
appear in the MHKIT must be integrated over 
ranges in which the Hankel functions do not have 
the sort of simple asymptotic form which would 
permit analytic simplification of our expressions. 

In such a case, the recourse is computation. The 
-> 
expressions for H were programmed onto a 7090 
computer. Typical results are shown for 2-ft waves 
at 20 kc/s in a "normal" sea 8 in figure 2 and for 25- 
ft waves in figure 3. The values of a were taken 
to be 7r/20, values typical of actual sea conditions. 
The figures give the attenuation (solid lines) and 
phase shift (dashed lines) of the magnetic field at a 
receiver located in a fixed horizontal plane, which 
plane is located 25 ft below the troughs of the waves. 
The data are plotted as functions of the height of 
the water directly above the position of the receiver 
in the horizontal plane. The portions of the curves 
drawn with thin solid lines represent points inter- 
polated from computed (solid and dashed) data. 

The reason for the existence of the vertical com- 

ponent of H, for propagation parallel to the wave 
crests, can be surmised by substituting H r from (7.2) 
for U in (6.4). This component is zero at any 
point directly under the crest or trough of a wave 
because of symmetry. It has a maximum value 
(for a fixed vertical coordinate) somewhere between 
the wave crest and the wave trough. 



8 Conductivity assumed 4 mho meter -*. 
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2-ft waves. 
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Figure 1. Coordinates for EM propagation over trochoidal 
cylinder. 
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Figure 3. Fields at a horizontal plane 25 ft under troughs of 
25-ft storm waves. 

8. Discussion of Results 

There are three points worth noting about the 
foregoing results : 

(a) With respect to the theory, we have assumed 
a simple trochoidal cylinder as the curve for the 
surface of the sea. However, the methods developed 
can be applied directly to any cylindric surface 
whose generating curve can be put in the form 

u=s—a 1 sin biS— a 2 sin b 2 s—a s sin 6 3 s— . . . 
v=a l cos b v s+a 2 cos b 2 s+ .... (8.1) 

In particular, a wave which is the sum of individual 
trochoidal wavelets has this form (provided the sea 
is not crossed). The static fields above such a sur- 
face can be obtained by a conformal transformation 
which is an extension of (2.2) : 

w=z+ja x exp jb x z+ja 2 exp jb 2 z + . . . (8.2) 

and so on, through all the formulas of section 2. In 
particular, the scale factor h is given by 
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h= Vl -22 ate-"*" cos b t s+^Za m a k e-^+V- cos (b k -b m )s. 

m, k 



If all the quantities a k are much less than one (as is 
usually the case for waves in the ocean), then this 
scale factor is approximately 



h=l- 



-XX^ & * n cos bjjs 

k 



(8.3) 



so that, roughly speaking, the superposition of waves 
results in a superposition of effects on the fields out- 
side the water. With respect to fields inside the 
water, the effects of superposing several waves are 
not clear a priori, since the boundary surface itself 
is altered. 

It is also worthwhile to remark that the MHIvIT 
for the underwater fields was obtained under the 
assumption that the external fields do not vary 

with the ^-coordinate; clearly, except for 0=90° 

-> 
(propagation perpendicular to the wave crests), H 
does have a ^-variation. However, the underwater 

propagation constant k is so large as compared with 

k sin 6 above the water that we are justified in ignor- 
ing the effects of a ^-variation in the integrand in 
arriving at (5.6) from (5.5). 9 

(b) With respect to the calculations, figures 2 
and 3 represent two extreme situations. In the 
former, the skin depth is large compared to the 
wave height and the depths of observation are 
comparable with the radius of curvature of the 
waves. Accordingly, the fields at depth show sub- 
stantial averaging-out of the variations. In the 
case of the 25-ft waves, the skin depth is less than 
the wave height and the depths of observation are 
small as compared with the radius of curvature of 
the waves. The fields in this case seem to behave 
substantially as if a flat sea were moving up and 
down over the point of observation. 

(c) With respect to possible experimental con- 
firmation of the theory, there are three possibilities. 
One can attempt to verify the predicted difference 

in the dependence of H on the configuration of the 
waves, either above or below the water surface. In 
figure 2, the phases of vertical and horizontal com- 

ponents of H depend differently on the amount of 
water above the observation point, and therefore 
vary differently with time; one could attempt to 
verify the elliptization thus introduced by the bound- 
ary. Finally, there is the difference in the phase 
behavior of the fields predicted at modest depths 
under relatively calm and strong sea conditions, as 
in figures 2 and 3. All three experiments are made 
uncertain by the fact that the real sea is not a simple 
trochoidal cylinder. The last two are also difficult 



9 T he effect o f assuming an e' k * dependence of the fields amounts to replacing \ 
by V(0H-(fr) 2 in (5.6). Since k«k, this change is small indeed. 



in that the effect to be measured is small . . . only 
a few degrees of elliptization or of phase change or 
a decibel or two difference in field strength. 

Experimental observations of VLF fields do, how- 
ever, confirm both the gross dependence of attenua- 
tion and phase on instantaneous wave height as 
predicted by the theory and also confirm the depend- 
ence of the local scale factor effect on the direction 
of atmospheric propagation with respect to the 
running of the sea. 10 The details of these experi- 
ments are described in the next section. 

9. Experimental Results 

The two conclusions of the theory that we put to 
experimental test are as follows: 

(a) For sea waves running from 5 ft swell to 
20 to 30 ft storm waves, the subsurface field is the 
same as if the sea were flat and the antenna were 
moving up and down, to within a few degrees in 
phase and a decibel or two in attenuation. 

(b) The E-tield in the air is locally enhanced over 
the crests of the waves and diminished over the 
troughs, in accordance with the local curvature of 
the surface. Whether this effect is passed on to the 
//-field and (through continuity of H at the bound- 
ary) to the subsurface fields depends on the angle 
between the direction of propagation and the direc- 
tion of the wave crests. 

To test the validity of conclusion (a) we made 
arrangements to tow a submerged float about 30 ft 
below the troughs of 20 to 30 ft storm seas in the 
North Atlantic, and to record the instantaneous 
variations in the magnitude (in decibels) and phase 
of VLF transmissions received on antennas in the 
float. 

Some of the results are shown in the strip chart 
recording of figure 4. The upper two traces are the 
received field strengths of VLF station NSS at 22.3 
kc/s on each of two underwater loops having orthog- 
onal directivity patterns. The next two traces are 
the corresponding variations in the received phase 
measured with a recording phasemeter. The bottom 
trace is a record of the variation in height of water 
above the towed float, as measured by a pressure 
transducer. Except that the pressure trace varies 
in the opposite sense, there is a clear excellent 
correlation between these curves, instant by instant 
in time. 

To properly observe the scale factor effect, con- 
clusion (b), requires better control over experimental 
conditions than is available in the open sea. To this 
end we made use of a field site at an old quarry wharf 
at Folly Cove in Rockport, Mass. 



10 The effect to be measured is the difference between /( and unity. From 
(2.7) the ratio of the maximum to minimum value of h is (l-fa)/(l— a). In a sea 
driven by strong winds a can rise from its normal value of 0.05 to 0.15 to the nearly 
breaking- waves value of 0.3 to 0.4. 
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Figure 4. Magnetic field and pressure measurements — open 

sea. 
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Figure 5. Folly Cove instrumentation. 

Part of the test equipment is in an underwater 
float in 70 ft of water 100 yards off shore. This 
includes a pair of loop antennas, preamplifiers, 
pressure transducers to measure depth, and telemetry 
as shown in figure 5. An anchor cable passes through 
a pulley to a winch on shore, so that the assembly 
can be raised and lowered. Other cables and ropes 
orient the loops and carry signals. On shore there 
are instruments for simultaneously receiving and 
recording two stations at once; they measure signal 
amplitudes in decibels and signal phase. A six- 



channel chart recorder on which to record and 
compare these measurements along with the cor- 
responding variations in head of water above the 
antenna loops is also used. The phasemeters oper- 
ate by comparing the received phase of standard 
transmissions such as those of NBA with that of 
a local precision frequency synthesizer. Since the 
transmissions of NBA and GBR are often keyed, we 
key the averaging circuits of the phase meter, so as 
not to measure the phase of the intervening noise. 
It has been found that besides those stations set up 
for standard frequency transmissions, many other 
VLF stations exhibit sufficient carrier stability to 
permit us to measure the phase of the received signal 
with our equipment. Thus, we have had greater 
flexibility in choosing signals than we might have 
originally expected. 

Since we measure depth with pressure transducers, 
we must take into account the fact that the pressure 
variation p due to the overhead passage of a given 
wave diminishes exponentially with the average 
depth d of the transducer. This and the hydro- 
dynamic relationship between the wavelength X and 
period T of gravity waves in open ocean are given by 

p=p e~ 2Td/ * 

A=1.56jT 2 meters, seconds 

p=p e~ id/T2 

in which d is depth of observation and p is the 
magnitude of the pressure variation as it would be 
at the surface. The last formula is a combination 
of the first two. The waves we see at Folly Cove 
do obey these laws. Apparently the presence of 
an ocean bottom at depth D (150 to 250 ft) for 
many miles simply eliminates by friction at the 
bottom those waves for which X is greater than the 
depth. The waves that remain are true gravity 
waves with sufficiently short wavelength as to have 
little interaction with the bottom. 

What the theory predicts is this : if the atmospheric 
propagation is across the sea wave crests, the under- 
water field will vary in attenuation according to the 
flat sea model as the waves go up and down. But 
if the atmospheric propagation is parallel to the 
sea wave crests, in addition to the flat sea effect 
there will be an enhancement of the field at the 
crests and diminution at the troughs of the waves. 
This effect operates in the opposite sense to the 
attenuation with depth, so as to reduce the total 
variation in attenuation. The phase shift, however, 
is not at all affected. If the direction of atmospheric 
propagation makes some general angle 6 with the 
direction of the wave crests, the theory predicts 
total underwater H field can be written as the sum 
of two orthogonal components, H x and i? n one of 
which varies with surface shape and one of which 
does not. 

In order to make any observation of this effect, it 
must be of the order of 4 to 6 dB, which in turn 
requires almost-breaking waves (wave height to 
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wavelength ratios of the order of 0.1 or greater) 
of a kind that exist only in a heavy storm. 

There was one storm at Folly Cove that produced 
almost ideal conditions for observing the effect, on 
a < [ayon which a NWgale blewfor 1 8 hours. The results 
are illustrated in figure 6. This is a strip chart 
recording in which the upper pair of traces are the 
received amplitudes (decibel scale) of signals trans- 
mitted from distant stations operating at 18.6 kc/s 
and 22.3 kc/s. The third and fourth traces are the 
corresponding received phases. The bottom trace 
is pressure (this time recorded in the same sense as 
the signal variations). By good luck, the 22.3 kc/s 
signal was propagating at about 45° with respect to 
the waves; moreover, the antenna loop directivity 
patterns were such that loop No. 1 saw mainly 
the component of H that is unaffected by surface 
shape and loop No. 2 saw mainly that component 
of H which is affected by surface shape. One would 
expect that the variation of amplitude introduced 
in the 22.3 kc/s signal in loop No. 2 by the surface 
shape would almost completely cancel out the 6-9 
dB variation of amplitude due to the varying height 
of water above the antenna (4 to 6 ft waves). One 
would also expect the 22.3 kc/s signal into loop 
No. 1 to show the full range of amplitude variation. 



SIGNAL AT 18.6 kHz 



SIGNAL AT 22.3 kHz 



BOTH SIGNALS 

^ FROM SAME - »[ ■ « 

LOOP (No. 2} 



1.3 kHz SIGNAL NOW 
FROM LOOP No.1 



PHASE AT 18.6 kHz 






BOTH SIGNALS 
— FROM SAME- 
LOOP (No. 2} 



PHASE AT 22.3 kHz 
! 



22.3 kHz SIGNAL NOW 
FROM LOOP No.i 



PRESSURE UNDER WAVES 



Figure 6. Magnetic field and pressure measurements — Folly 

Cove. 



Tn the left hand part of the records in figure 6, 
both receivers were operating from loop No. 2; 
in accordance with theoretical expectation, the 
amplitude of the 22.3 kc/s signal shows little or no 
variation. Midway through the record, the 22.3 
kc/s receiver was switched to loop No. 1 ; as pre- 
dicted, the full 6 to 9 dB range of signal variation 
is apparent. Note that the phase variations are 
the same in both cases. 

As a matter of fact, the influence of the scale 
factor effect on the field observed in a directive 
underwater antenna (e.g., a loop) varies in a compli- 
cated way with the angle between the direction of 
atmospheric propagation and the wave crests, the 
orientation of the antenna, and its pattern. It is 
significant that by taking these factors into account 
it was possible to calculate the approximate direction 
of the two stations being received in figure 6 from 
the known orientation of the loops and the changes 
in average signal strength and phase as the receivers 
were switched between loops. These calculated 
directions agreed within ±5° of the directions 
derived from the known locations of the trans- 
mitters. It is even more significant that taking 
account, in addition, the observed changes in the 
range of amplitude variations as the signals are 
switched between loops, it is possible to calculate 
the direction of the running of the sea. r I nis calcu- 
lation agrees (within 30°) with what was observed 
as closely as it was possible to estimate that direction 
by eye and hand compass. 

10. Appendix A 

We wish, here, to find a solution of the scalar 
Helmholtz equation in two dimensions 



VV(«,fO+JfcV(«,n)=0 



(Al) 



in order to obtain a scalar function <p from which to 

-> 
derive an E'-field 



2?= curl cpe 2 



(A2) 



where E is the electric field associated with an 
electromagnetic wave which propagates in a direc- 
tion normal to the sea wave crests. The tangential 

-> 
component of E must vanish at the surface n=0. 
Since, in the limit k=0, the fields must approach 
the static field pattern, we begin by taking 



<p=\p -- 



jk 



(A3) 



If we place this ^ in the Helmholtz equation (4.1) 
or (4.9), we obtain 

vVo+&¥o= J° V (l -£) e~ m (A4) 
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where since we consider propagation perpendicular 
to the wave crests, sin 0=1. Thus in the limit fc^O, 

VVo=0 

satisfying Laplace's equation, and 



M=.Z?=curl fae z =-T E e n . 



(A5) 



This is precisely the static field that results from 
conformal transformation of the field between 

parallel plates. It is unusual only in expressing E 
(in a charge-free region) as the curl of a vector 
potential rather than as the gradient of a scalar 
potential. When k^O, fa satisfies the scalar wave 
equation only if h 2 = l; in this case, the coordinate 
system would be rectangular and fa would represent 
a plane wave traveling in the ^-direction between 
parallel-plate bounding surfaces. 

In the present problem, h 2 ^l, the right-hand side 
of (A4) is not identically zero, so that fa is not an 
exact solution to the wave equation in trochoidal 
coordinates. Suppose, now, that we can write a 
second approximation to the solution of the wave 
equation as 

fc=*,+^«-* (A6) 



in which Q x is a function of n and s for which 

V 2 &=-F(l-1/A 2 ) (A7) 

and in which Q x is a function whose maximum mag- 
nitude is of the order of magnitude P. Substituting 
fa in the Helmholtz equation, we obtain a result 
which has the form 



where 



Wi+*Vi=§ PM s)e~ jks (A8) 

jfc 

P x (n, *)=[2j* ^+F(l-1/A 2 )ft]. (A9) 



Pi is a function of the order of magnitude k z in- 
stead of one of the order of magnitude k 2 as in the 
case of fa and (A4) . Once more, we can take a new 
trial solution \f/ 2 . 



fa=fa 



Q2E0 
jk 



- jks 



(A10) 



in which Q 2 is a function of n and s for which 



(All) 



and in which Q 2 is a function whose maximum 
magnitude is of the order of magnitude P! This in 
turn leads to an equation which has the form 



V 2 fa+k 2 fa=^P 2 (n, s)e~ jks 



(A12) 



in which P 2 is a function of the order of magnitude 
& 4 instead of F as for the case of fa. Thus, we are 
led to a sequence of functions 

fa, fa, fa, . . ., fa 

which satisfy the Helmholtz equation in the limit as 
n is taken sufficiently large, 

UmV 2 fa+^fa=0 &<1 



provided that k is less than one. At the same time, 
we have found a solution to the scalar Helmholtz 
equation 



7l->oo 



jk 



is of 



in which the maximum amplitudes of the Q m 
the order of magnitude k m ~ l times that of Q lm 

The original trial solution fa satisfied the boundary 
-> -> 
condition that M=E be normal to the water surface, 
but did not satisfy the Helmholtz equation. At 
each step of the approximation, it is necessary to 
pick the #'s so that \f^\, fa, . . . continue to satisfy 
the boundary condition. It is beyond the scope of 
this paper to show that in the differential equations 
for the Q's there is sufficient flexibility to carry 
through the program just outlined. 11 

Now, let us carry out the program for the case at 

hand. Our fa is given by (A3). When substituted 

-> 
in (A2), ^ leads to an I^-field normal to the surface 
n=0, but does not satisfy the Helmholtz equation 
(Al). Using the expression for V 2 in (3.1) and that 
for h in (2.7), we find that in trochoidal coordinates 
(A7) for Qi becomes 



a 2 «i , v 2 Qa_ 



ds 2 i>n 



-2k 2 ae~ 



cos s 



-We-z* 



A particular solution of (A13) is 

a p a 2 

Qi=—k 2 ane n cos s — -r— e 2n . 



(A13) 



(A14) 



To this solution, we can add any number of terms 
which are solutions of the homogeneous equation 
associated with (A13) 



Qi=A r e rn cos rs 



(A15) 



in which r is any real number and A r any complex 
number. In general, terms from the set { Q x } can 



11 Indeed, in order to satisfy the boundary conditions exactly at each step, one 
must generally add to Q m a term which is an approximate solution of the Helm- 
holtz equation rather than being a solution to a differential equation like (A7). 
Such terms have sufficiently small coefficients so that they do not upset the 

convergence of the sequence of functions, ^1, fa They amount to slight 

modifications of fa, having relative magnitudes less than k m , based on the re- 
quirements of the mth step. 
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be added to Qi for the purpose of satisfying the 
boundary that E be normal to the surface n=0. 
We remark that both terms of $1 lead to tangential 
components of !£ at n=0. In order to eliminate 

the tangential Z£ at n=0 due to the first term on the 
right of (A14), we add to (A14) a term from the 

set { Qi } to obtain for Q x 



Q l =—k 2 a(n-{-l)e w coss- 



k 2 a 2 



(A16) 



The tangential components of E due to the second 
term of (A14) cannot be counteracted by terms 
from the set (A15) with r=0, because such terms 
do not meet the reasonable requirement that their 
effect tend to zero as n gets very large. To com- 
plete the satisfaction of the boundary condition, 
we add to \pi a further term pi. 



ka 2 



<Pi=- 



jk 



jks. 



(A17) 



Thus, the modified second trial solution to the 
Helmholtz equation (Al) becomes 



= l—k i a(n-\-l)e~ n cos s- 



k 2 a 2 



ka 2 _,„"|£ 



-jks. 



(A18) 



This is the function \f/ given in the text in (4.11). 
If we use the Helmholtz operator on \p, we obtain 



v-v+fcv 



-k 3 [ 2 ja (n+l) e~ n sin s+^ e~ kn ~] 
+ terms of higher order in k. 



(A19) 



Thus, \p is a solution of the Helmholtz equation to 

the third order in small quantities, k. 
-> 
To derive the Z?-field from (Al8), it is convenient 

to use the vector identity which states that if u 

— > 
is a scalar function and F is a vector function, 

respectively, of coordinates, then 



V X (uF) = (Vu) X F+ uV X F. 



(A20) 



Thus, for (A2) we have 
E=cur\ \l/e z =VXtye Zd 






)Xe z 



:> 



(A21) 



We can thus use (A21) to compute the /t-field due 
to the scalar solution determined in (Al8). The 
operations are straightforward and lead to the result 



jkh 
jkh 



[^ ( «-- 



e kn )-\-k 2 ane n cos s \e 



} 



* , \jk-k i (a[n+l]e- n sin s)-jk*(jt 



— — e kn -\-a[n + l]e "coss 



(A22) 



This is the same expression as that given in (4.13) 

of the text, except for the cancellation of the factor 

-> 
jk. The tangential component of E vanishes at 
n=0; the tangential terms in e~ n and e~ 2n come from 
corresponding terms in Q { ; while the term in e~ kn 
comes from <p. 

11. Appendix B. Reduction of the H—K 
Integral 

We wish to reduce the Helmholtz-Kirchoff integral 
which appears in (5.5) by integrating with respect 

to z. 



u ^=&[ u ^-"°Ut) 






] 



dS. (5.5) 



We let a plane z=z intersect the cylindric surface in 
a curve C, which has the same shape regardless of 
the value of the constant z . Let us denote the 
differential element of arc length on C by ds; then 
dS=dsdz and we can rewrite (5.5) as 



_j_ r dU(u,v) r r 

4w Jc cto |_J - 






] 



e~ jbz dz Ids. (Bl) 



Further, if the point P lies in the plane z=z , we 
can set £=-y/(z— z ) 2 +r 2 , where r is the length of the 
projection on the plane z=0 of the line from P to a 
point on the curve C. 
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So we have 



;/cV(2-z ) 2 +»- 2 



If we make a change of variables u=z— z , we get 



/. 



#V« 2 +r 



V(3-3 ) 2 +r 2 
We note that 



e- ibz dz=e- Jbz ° . e- jb »dcc. 



(B3) 



jk-yj^+r 



Vco 2 +r 2 
is an even function. Hence 



co 7'/ c y a ,2-|_ r 2 



J"; 

«y — co -y 



/,.,2_)_ r 2 



<?-^ w d. 



H". 



Vco 2 +r 2 



cos 6coc?w. (B4) 



From well-known tables, 12 we have 

;/cVco2+r2 



J -co A 



/,., 2 -fr 2 



cos b^dcc=J7rm i) [r(k 2 -b 2 y /2 ]. (B5) 



Therefore, setting (lc 2 —b 2 ) l,2 =ki i becomes 
U(P)=j e -^j c \U{u, v) ± fly for) 

We will note also that if we choose P and C in the 
same plane perpendicular to the cylinder axis, the 
factor e~ jbz ° disappears and we have 



U(P)=l j c [u(u,v)^H^(t 



r) 



-H^ikr) 



dn 



1 ds. (B6) 



Now if r is the radius vector from a point p on 



12 "Tables of Integral Transforms," Bateman Manuscript Project (McGraw- 
Hill Book Co., Inc., New York, N.Y., 1954), Formulas 1.7 (34) and 1.7 (30) on p. 26 

and 1.5 (27) on p. 17. We need also Ho<» (z) =Jo(z)-\-jYo(z) and Ko(z) =^ W» 

( ze -J j which may be found in "Higher Transcendental Functions," Bateman 

Manuscript Project (McGraw-Hill Book Co., Inc., New York, N.Y., 1953), Vol- 
ume II, as Formulas 7.2.1 (5) and 7.2.2 (15) respectively (pp. 4 and 5). 



C to P and n is the inward unit normal to C at p t 
then 



dn r 



(B7) 



U(P)= 3 -jJjcUW»(kr) r -^-HP(Jcr) g] ds . (B8 ) 

This result is the Modified Helmholtz-Kirchoff 
Integral Theorem (MHKIT) for cylindric geometry. 
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